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CHAPTER I 


INTRODUCTION 

Chapter One presents the background and overview of this study. A 
brief review of the lower atmosphere and its electrical parameter, the 
complex microwave refractlvity, will be given. 

A. Overview 

As the use of the electromagnetic spectrum has expanded beyond ten 
GHz, a number of theoretical and experimental studies on the effect of 
the atmosphere on propagating waves have been performed, and the atmos- 
phere has been found to have a strong influence in this portion of the 
spectrum. In this high frequency region, the refractive index of the at- 
mosphere can no longer be regarded as a frequency independent quantity 
in discussions of its influence on radio wave propagation. Microwave si g 
nals travelling along atmospheric paths undergo attenuation and refrac- 
tion due to absorption by atmospheric gases and variations in the refrac 
tive index, respectively. After Van Vlecks's pioneering work on micro- 
wave absorption theory (1); several models were introduced to develop 
the microwave refractivity of the clear sky. Based on these works, the 
microwave refractivity is highly dependent on temperature, pressure, and 
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and water vapor content which, in turn, are random functions of time and 
space. 

The design of reliable earth-space and terrestrial communication 
links requires estimates of the deterministic and random changes of 
these atmospheric effects on radio wave propagation. 

Using the extensive meteorological data base available from the Na- 
tional Weather Service, the diurnal and seasonal variability of the sur- 
face microwave refractlvity has been determined at Columbus, Ohio. The 
results of this effort make it possible to estimate the variability of 
the absorption due to atmospheric gases and the refraction effects for 
a standard atmosphere. 

In the following, a standard exponentially tapered atmosphere is 
assumed and the effects of turbulence are Ignored. 

B. Atmosphere 

The atmosphere that surrounds the earth is a mechanical mixture of 
many gases. Although the atmosphere is composed of a number of gases, 
five of them - nitrogen, oxygen, argon, carbon dioxide, and water vapor - 
make up 99.997 percent of it by volume below 90 km. 

A complete description of the physical state of the atmosphere at 
any moment requires the measurement of dozens of quantities. The analy- 
sis of weather is usually centered on several variable characteristics 
of the atmosphere near the earth's surface: temperature, pressure, humid- 
ity, cloud, precipitation, and winds. Among these characteristics temp- 
erature, pressure, and humidity are major factors in the determination 
of the refractive index. Figure 1 shows the vertical distributions of 
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Figure 1. Characteristics of U.S. Standard Atmosphere, 1962, compared 
to those of ARDC Model Atmospheres, 1956 and 1959. (The 1956 
Model agrees with the 1954 "U.S. Extension to the ICAO Stand- 
ard Atmosphere".) The names shown ^or atmospheric shells 
and boundaries are those adopted by the World Meteorological 
Organization. 
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temperature, pressure, and density for the U.S. Standard Atmosphere. 

1 . Temperature 

Temperature Is a measure of kinetic energy of the molecules pre- 
sent In the air. Temperature varies widely from place to place over 
the earth's surface as well as through time at a given place. There 
are two major periodic changes in temperature: the annual and the di- 
urnal . 

As shown in Figure 1, the temperature normally decreases with 
height in the troposphere. An approximate formula which describes the 
temperature-height dependence in the troposphere is (2): 

T = 229 - 6.55z [°k] ; 0 < h < 11.6 km 

O-l.l) 

= 233 [°k] ; 11.6 < h < 20 km 

Where T = temperature at z 

z = height above mean sea level [km] 

2. Pressure 

The pressure exerted at any level is due almost entirely to the 
weight of the air pressing down from above; this force results from 
gravitational attraction. As would be expected, the pressure changes 
rapidly in the vertical direction. The rate of its change is greatest 
near the surface and decreases at a decreasing rate with height. The 
pressure-height dependence can be approximated by an exponential law (3) 

P = P 0 exp [- z/H ] [mbar] (1-1.2) 

Where, 

P = total pressure at altitude, z, 

P 0 = sea level pressure averaged over the earth 
y 1,013 [mbar] 

'V 
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H = scale height 
n, 7 [km] 

z = height above mean sea level [km] 

Variations of pressure in the horizontal are much smaller than 
those in the vertical. Near sea level, the change of pressure with dis- 
tance rarely exceeds 3 mb per 100 km and usually much less than half of 
this rate (4). Only In the presence of severe storms or exceptionally 
well developed masses of cold air does the pressure vary more than three 
percent from the average. Even though the variation in the horizontal 
Is relatively small, this small change in pressure may change the direc- 
tion and velocity of the wind and can subsequently lead to temperature 
and moisture changes. 

Pressure also changes with time at a single place. Some of these 
changes are of an irregular nature, caused by occasional invasions of 
air having a different mean density. But there is also a quite regular 
diurnal oscillation of the pressure. 

3. Water Vapor 

The concentration of water vapor in the atmosphere varies from 
practically zero to as much as four percent by weight. The vertical 
profile of the averaged density of water vapor, pw, can be approximately 
described by an exp nential law (2) in a humid atmosphere, 

(^ = 20 exp [-0.1 2z] [gm/m 3 ] (1-1.3) 

while, for a standard atmosphere, 

p u = 5 - j [gm/m 3 ] (1-1.4) 

The parameter which is in common use for the representation of the 
water vapor content of the atmosphere is the relative humidity. The 
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definition of relative humidity, v, is: 

v = j- X 100 [i] (1-M; 

c s 

where e is the partial pressure due to water vapor and e $ is the satura- 
tion pressure. Starting with the latent heat of vaporization, an ex- 
pression for the saturation pressure may be derived with the following 
result: 

e s = 6.11 exp [(6,830 - 5.14 T) {^y - y)] [mbar] (1-1.6) 

where T is temperature [°k]. 

C. Refractive' ty 

The mathematical formulation for the electric field E of a uniform 
plane wave propagating in a 1-direction can be written: 

E = E 0 exp [-jklj (1-2.1) 

where k is the complex wave number. From Jfexwell's equations, 

k = u = ^-y ~T r (1-2.2) 

where t r is the complex relative permittivity. The complex index of re- 
fraction of the atmosphere is then giver, by: 



where the prime and double prime denote real and imaginary parts, re- 
spectively. Usually the complex microwave refractivity, N, is used for 
convenience in practical acrlica.ions, 

N = (n - 1 ) 10 6 = N'- jN M [ppm] (1-2.4) 

There is another propagation parameter in comnon use: that is the propa- 
gation corstant: 

Y = jk = a £ + j S (1-2.5) 
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where, s E is the specific attenuation of the electric field and g Is the 
phase constant. F-orn equations (1-2.2), (1-2.4) and (1-2.5), 


<* E * 0.02096 fN" [neper/km] 0-2.6) 

S = 0.02096 f(10 6 + N' ) [rad/krc] (1-2.7) 

wt«re f is in [GHzj. The specific power attenuation, a p , can be defin- 
ed as: 

= 10 loq j a J-*10 lo ?H (exp[0. 02096 x 2fN"l} = 

P ' P ° (1-2.8) 


0.182 fN” [dB/kro] 
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CHAPTER II 


METEOROLOGICAL DATA 

The weather data used in this report was measured by the National 
Weather Service at Cclunt>us, Ohio, from 1948 to 1976. Temperature, 
pressure, relative humidity, wind, and several other meteorological 
states were observed every hour during the 1949-1958 period, and every 
three hours during the 1965-1976 period. 

The data are stored on two nine-track tapes. The record format 
consists of fixed blocks each composed of four equal length logical re- 
cords. Each logical tape record contains six observations and is 495 
bytes long. These records consist of fifteen bytes of identification 
followed by six observations of eighty bytes each. Records always be- 
gin with the Local Standard Time hour of 00LST, 06LST, 12LST, or 18LST. 
Space is always etained on tape for twenty- four observations/day. When 
no observation is available, the hour is indicated (two bytes) and all 
other fields are coded blank. 

For more detailed information about the format and coding, the 
reader is referred to Appendix A. 
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CHAPTER III 

STATISTICS OF TEMPERATURE, PRESSURE, AND WATER VAPOR CONTENT 

Columbus, Ohio (10), is located in the area of changeable weather. 

Air masses from central and northwest Canada frequently invade this re- 
gion. The tropical Gulf air masses often reach central Ohio during the sum- 
mer and, to a much less extent, in the fall and winter. There are also 
occasional weather changes brought about by cool outbreaks from the Hud- 
son Bay region of Canada, especially during the spring months. At in- 
frequent intervals the general circulation will bring showers or snows 
to Columbus from the Atlantic. Although Columbus does not nave a "wet" 
or "dry" season as such, the month of October has a higher frequency of 
light rainfall than any other month and comes closest to providing a nor- 
mal dry period. 

In this chapter, the statistics of temperature, pressure, and water 
vapor content are obtained for two ten-year periods: 1949 through 1958, 
and 1967 through 1976. Note that the weather data were observed every 
hour during the 1949-1958 period, and every three hours during the 1967- 
1976 period. 

For each weather variable, the probability density, the average di- 
urnal variation on a monthly basis and the daily maximum, average, mini- 
mum, and extremes by months are clotted. The maxima and minima are the 
averages of each year's highest and lowest values. The extremes are the 
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extreme values observed during the period specified. 


A. Temperature 

Since the transient periods between summer and winter are relative- 
ly short and the two seasons are distinctive, the probability density 
functions of temperature have a bimodal frequency curve with two maxima 
at about 31°F and 70°F (see Figures 2 and 3). The diurnal variations of 
temperature duri-.g summer are greater than those of winter (see Figures 
4 and 5). But the fluctuation during winter is a little greater than 
that during summer (see Figures 6 and 7}. In all cases, the temperature 
reaches its minimum value near local sunrise and its maximum at about 
3:00 p.m. local time. 

B. Pressure 

The probability density function has the shape of Gaussian distrib- 
ution (Figures 8 and 9), and has a peak at 987.14 [rt>] in both periods. 
This value is slightly higher than the value of 978 [nfc] predicted by 
Equation (1-1.2) for the observing station height of 247 [m] above mean 
sea level. 

For the 1949-1958 period, the average is 987.2 [mb], the highest is 
1015.62 [mb], and the lowest is 949.24 [mL]. For the period 1967-1976, 
the average is 987.5 [mb], the highest is 1013.25 [mb], and the lowest 
is 953.65 [mb]. 

Figures 10 and 11 show that the diurnal variation of pressure is 
quite small. But it is interesting to note that there is a quite regu- 
lar diurnal oscillation of the pressure that causes, on the average, two 
maxima (at about 10:00 a.m. and 10:00 p.m.) and two minima (at about 
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Figure 3. Probability density function of temperature for the period 
1967-1976; U = 2[°F]. 
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Figure 4. Average diurnal variation of temperature for the period 
1949-1958. 
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Fiqure 5. erage diurnal variation of temperature for the period 
. '7-1976. 
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Figure 6. Daily maximum, average, minimum, and extreme temperature 
by months for the period 1949-1958. 
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Figure 7. Daily maximum, average, minimum, and extreme temperature 
by months for the period 1967-1976. 
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Figure 8, Probability density function of pressure for the period 
1949-1958; Ap= 1.43 [mbar]. 


1 




PROBRBIL I TY DENSITY 
L 0.00 0.02 0.01 0.06 0.08 


PERIOD : 1 /I /67-12/3i/76 



PRESSURE (MBRR) 


Figure 9. Probability density function of pressure for the period 
1967-1976; AP= 1.43 [mbar]. 
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Figure 11. Average diu.-nal variation of pressure for the period 
1967 - 1976 . 


4:00 a. m. and 4:00 p.m.). The dally heating and cooling cycle appears 
to be the dominant cause of the pressure variations. And, the pressure 
during the spring season tends to be lower than that during the autumn 
season. 

Even though the average pressure does not vary much seasonally, 
the degree of fluctuation of pressure is dependent on the season (see 
Figures 12 and 13). The fluctuation of pressure during winter is much 
greater than during 'umner. 

C. Water Vauor Content 

The probability density functions of the relative humidity (Figures 
14 and 15) show a strong bias toward values of high relative humidity. 
The diurnal variation is small during the winter season and large during 
the summer season (Figures 16 ard 17), but the monthly average is the 
highest in January and December (Figures 18 and 19). 

Since the water vaoor pressure is dependent on temperature (Equa- 
tion 1-1.3 and 1-1.4), the shape of its density function (Figures 20 and 
21) is quite different from that of relative humidity. The water vapor 
pressure doesn't vary diumally, but the monthly average shows a strong 
seasonal dependence (Figures 22, 23, 24, 25). The average in July is 
about 20 [mb], and the average in January is about 5 [mb]. 
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Figure 13. Daily maximum, average, minimum, and extreme pressure 
by months for the period 1967-1976. 
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Figure 14. Probability density function of relative humidity for 
the Defied 1949-195?; i v * ![?]. 
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. Average diurnal variation of relative humidity for the 
period 1967-1976. 
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Figure 19. Daily maximum, average, minimum, and extreme relative 
humidity for the period 1967-1976. 
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Figure 21. Probability density function of H^O pressure for the 
period 1967-1976; A e = 0.71 [mbar]. 
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Figure 22. Average diurnal variation of H ? 0 pressure for the period 
1949-1958. 


32 


H20 PRESSURE (MBRR) 



HOUR 


"uure 23. Average diurnal variation of H^O pres.ure for the period 
1967-1976. 
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Figure 24. Daily maximum, average, minimum. and extreme F^O 
pressure by months tor the period 194S-1958. 
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CHAPTER IV 


ATMOSPHERIC REFRACTIVITY VARIATION 

Atmospheric gases having electric or magnetic dipole moments cause 
absorption and dispersion of electromagnetic waves. In atmospheric pro- 
pagation, t^e atmospheric effects are represented by the complex refract- 
ivity, N, in parts per million [ppm]. Since the molecular polarization 
depends upon temperature, pressure, and gas composition, the refractiv- 
ity varies with the atmospheric state. 

In this chapter we will calculate the clear sky refractivity fer 
the meteorological conditions presented in the preceding chapter and its 
probability distribution and variation for Columbus, Ohio. 

A. Calculation 

Over the 10-300 GHz frequency range, the spectra of the atmospheric 
gas molecules lead to a complicated expression for N with several line 
parameters. The Absorption Line Parameter Compilation (ALPC) (5), which 
is the most complete account of atmospheric absorption lines, lists more 
than 10,000 l^nes for the seven molecules: H^O, CO^, 0^, N 9 0. CO, CH^, 
and 0 o . Of these, six water vapor and thirty-six oxvqen spectral lines 
dominate the electromagnetic characteristics of clear air in the EHF 
range. 

Based on Omoura and Hodge’s work (6), the complex microwave refract- 
ivity, N, can be expressed as the sum of the refractivities due to each 
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individual absorption line. 


n - rr^ = zn' - jirfq = n<, + r(N* - n c1 ) - « n c ♦ 

i=l 1=1 1=1 1 1 


i{S 1 F; - N..) - jtS ^ 1 


(4-1.1) 


where 

N 0 = rN 0 , is the frequency Independent real part of the refracvlv- 
1 1 

Ity of air [ppm], 

N 0 .j - the low frequency limit of the 1-th line refractivity in 
[ppm], 

Sj = the i-th line strength In [KHz], 

F^' = the dispersion line shape factor of the 1-th line in 

[GHz”^]. and 

F^" = the extinction line shape factor of the i-th line in 

[GHz' 1 ]. 

As will be shown in Section ( A-l ), the low frequency tail can be ap- 


proximated as: 

S i 

No, = lim S F! = S. 11m F ! = 2 *- 
1 f+° 1 1 1 f^ e 1 


(4-1-2) 


where is the resonance frequency of i-th line. Therefore, if we 
neglect the effect of absorption lines other than the six water vapor 
and thirty-six oxygen lines, the clear air refractivity in the EHF range 
can be obtained by 

( 42 S. \ 6 36 

N °\:I, * if, (s i F ;>wt * (s i F i>ox < , - 1 - 3 > 
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( 4 - 1 . 4 ) 



< S i F i> wt 


36 

♦ r (s.Fy) 

i=l 1 ‘ 


ox 


where "wt" stands for water vapor and "ox" stands for oxygen. 

The oxygen spectral lines and the water vapor spectral lines with 
their line parameters are listed In Tables 1 and 2. From Equation 
(1-2.8) the power attenuation will be given by 

« = 0.1P2f ^ (S t Fp. t * 
where f Is in [GHz]. 


1 . Frequency Independent Refractlvity 

The frequency independent refractivity can be obtained by summing 

the low frequency limits of all of absorption line refractivities, i.e. 

N 0 = I N . i = all absorption lines (4-1.6) 

i 01 

But this is not possible because it requires a full knowledge of all of 
the absorption lines. Based on experimental work. Bean and Dutton (7) 
presented the following empirical equation for the frequency independent 
refractivity of moist air: 

N 0 = 77.6 ^ + 72 j + 3.75 x 10 5 ^2 (4-1.7) 

where 

D d = dry air pressure [mb] 
e = partial pressure of water vapor [mb] 

T = temperature [°k] 

Note that this frequency indeoendent refractivity is purely real. 

As shown in the next section, the imaginery part of the refractivity is 
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TABLE ] 

OXYGEN SPECTRAL LINES (5) 



* [■*-'] 

S 3 [kHz/*b] 

d 

> 0 [GHz/«b] 

1 

5.5 xl0~ 4 

2.44 xlO' 7 

8.69 

8.7 x 10“ 4 

1 

5.6 x 10' 4 

6.04 x 10 -7 

7.74 

8.9 x IQ' 4 

53.02117 

5.5 x 10" 4 

1.41 x 10" 6 

6.84 

9.2 x 10 -4 

52.54223 

5.7 x 10’ 4 

3.08 x 10' 6 


9.4 x 1 0 -4 | 

53.0668 

5.3 x 10 -4 

6.37 x 10’ 6 

— 

5.22 

9.7 x 10’ 4 ^ 

53.59572 

5.4 x 1C' 4 

1.24 x 10' 5 

4.48 

10.0 X 10- 4 i 

54.12997 

4.8 x 10' 4 

msm 1 

3.81 

10.2 x 10' 4 i 

54.67116 

4.8 x 10' 4 

HI 

3.19 

10.5 x 10‘ 4 ; 

55.22136 

4. 17x 1(T 4 

6.274x 10' 5 

2.62 

10.79x lO' 4 ! 

55.78380 

3. 75x 10‘ 4 

9.47lx 10‘ 5 

2.11 

ll.lOx 10~ 4 i 

1 

56.26478 

7. 74x IQ* 4 

5.453x 10' 5 

0.0138 

16.46x 10 _4! 

56.36339 

2.97x 10* 4 

1.335x 10‘ 4 

1.66 

11 44x 10' 4 | 

56.96818 

2. 12x 10' 4 

1.752x 10‘ 4 

1.255 

11. Six 10 4 f 

| 57.61249 

0. 94x 10‘ 4 

2. 126x 10’ 4 

0.910 

12.66x 10' 4 

58.32389 

-0.55x 10‘ 4 

2.369x 10- 4 

0.621 

12.66x 10“ 4 ; 

58.44660 

5.97x 10' 4 

2.0827xl0 -4 

0.0827 

14.49x 10~ 4 ; 

59.16422 

-2.44x 10' 4 

2.387x 10' 4 

0.386 

13.19x 10' 4 i 

59.59098 

3.44x 10" 4 

2.097x lO' 4 

0.207 

13.6 x IO- 4 ! 

60.30604 

-4.35x 10* 4 

2.109x 10* 4 

0.207 

13.82x 1C' 4 

! 60.43478 

i 

1.32x 10‘ 4 

2. 444x 10‘ 4 

0.386 

12.97x 10‘ 4 

61.15057 

-0. 36x 10‘ 4 

2.486x 10' 4 

0.621 

12.48x 10" 4 

61.30017 

-1.59x 10' 4 

2.281x 1C* 4 

0.910 

l?.07x lO' 4 

62.41122 

-2.66x 10' 4 

1 .S19xl0~ 4 

1.255 

1 1 . 7 1 x lO" 4 

62.48626 

-5.03x 10‘ 4 

1.507x 10' 4 

0.0827 

14.68x 10" 4 

62.99800 

-3. 34\ 10* 4 

1.492x 10 -4 

1.66 

11.39x 10' 4 

63.56852 

- 4 . 1 7 x 10‘ 4 

1.079x 10' 4 

2.11 

O 

00 

X 

0 

1 

64.12778 

-4.48x 10‘ 4 

7.281x 10' 5 

2.62 

10.78x 10‘ 4 

64.67892 

-5.1 x 10' 4 

4.601x 10 -5 

3.19 

10.5 x 10 -4 
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Table 1 (continued) 


6:;. 2240* 
65. 76474 

BS HEM 

RnVhU 

BSZES1 

3 

3.81 

4.48 


6i> . 302Qo 
66.83677 

-5.5 x lO’ 4 
-5.9 x 10 -4 

7.94 x 10’ 6 
3.91 x 10~ 6 

5.22 

6.00 

9.7 x 10' 4 
9.4 x 10' 4 

l 

67.36951 

67.90073 

-5.9 x 10' 4 
-5.8 x 10' 4 

1.81 x lO' 0 
7.95 x 13~ 7 

6.84 

7.74 

9.2 x 10~ 4 
8.9 x 10' 4 

68.4308 

118.76034 

-5.7 x lO’ 4 
-0.44x 10" 4 

3.28 x 10~ 7 
9. 341x 10* b 

8.69 

0.C138 

3.7 x 10‘ 4 
15.92x 10' 4 


Table 2 

WATER VAPOR MICROWAVE SPECTRAL LINES (5) 


BBI 

S 0 kH 2 /mbar 

— 

d 

> 0 GHz/mbar 

22.23508 

0.0112 

2.143 

28.1 x 10‘ 4 

63.052 

0.0018 

S. 75 

28.0 x iO' 4 

183.31009 

0.241 

0.653 

28.2 x 10' 4 

321.22564 

0.0044 

6.16 

22.0 x 10* 4 

325.15292 

0.159 

1.52 

29.0 x 10' 4 

380.19737 

1.240 

1.02 

28.5 x 10~ 4 
L 
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proportional to the operating frequency for low frequencies and Is gen- 
erally negligible for frequencies below 1 Wz. 


2. Line Shape Factors 

The frequency variation of the absorption and dispersion of each of 
the spectral lines Is represented by a line shape factor. Three classi- 
cal line shapes (the Lorentzlan, the Van Vleck-Welsskopf , and the Gross) 
are In common use In atmospheric refractlvlty studies, but the Gross 
line shape seems to fit the experimental results better. The three line 
shape factors and their tails are listed In Tables 3 and 4, respective- 
ly. 

At low altitudes (h < 35km), the effect of pressure broadening 
causes the absorption lines to overlap. Based on Rosenkranz's work, 
Llebe (5) used the following line shape factors for the calculation of 
clear air refractlvlty In the pressure broadening region (h < 35km). 


"<9 
- <fr) 


(to - f) + yl 

(fo - f) Z + y 2 


Y - ,(f, - f) I 


(f» - f) 2 + y 2 


(f 0 + f) + yI 


Y - -(t , - t . f ) l 

TTT2 


(4-1.8) 


(4-1.9) 


(fo + f f + r 

where f„ Is the resonant frequency of the line, f is the operating fre- 
quency, y Is the line width, and I is an interference parameter. The 
Interference parameter is negligible In the case of water vapor while It 
is a function of pressure and temperature for the oxygen molecule (5). 

2 

(4-1.10) 


I = a Pd 


(300) 


where 
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TABLE 3 


LINE SHAPE FACTORS (6) 


Type 

Dispersion F' 


Full Lorentzlan 
(F-L) 

fo-f , fo+f 

Y Y 

(fo-f) 2 -*y* (fo+f) Z +v 2 

(fo-f)^Y Z (fo + f) Z+ Y 2 

V*n Vleck-Welsskopf 
(VV-W) 

/V+Uf.+f) , v 2 +fo(fo-f)V IN 


\Jf 0 +f) 2 +Y 2 (f.-f)V j\ V 

(f.-f)Sy v T v 

Gross 

(6) 

l 

4f 2 Y 

[f l - f 2 ] 2 + 4fV* 





TABLE 4 


LINE SHAPE FACTOR TAILS (6) 


m 

f-TQ 

F' P 

f -► « 

F* F" 

1 

2 

77 

4 Y f 
f ^ 

1 O 


4yf 0 

~r~ 

vv-w 

2 

77 

f J 

* o 

O fo 

& 

G 

2 

fo 

A" 

» o 

-2 % 
r 

$ 
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a = constant 


Pd = dry air pressure [mb] 

T = temperature [°k] 

If y << f„ and I is small, the low and high frequency tails are 

F' % ; F 11 % as f ->■ 0 (4-1.11) 

f 0 f J 

I o ' o 


F' 




as f -*• ® 


(4-1.12) 


For comparison, we plot the dispersion and extinction line shape factors 
of the four lines in Figures 26 and 27 for equal tu 0.065. The in- 
terference parameter is set to 0.06. As shown in Figure 26 and from 
comparison of low and high frequency tails, liebe's dispersion line 

A 

shape factor seems to be shifted downward by 4— . It should also be no- 

* o 

ticed that Liebe extinction line shape factor, F", goes below zero if 
j- < I [see Equation (4-1.11)]. Several oxygen absorption lines satisfy 

* O 

this condition for usual atmospheric states. For these two reasons, the 
Gross line shape will be used in the following rather than Liebe's more 
recent result. 

The line width and line strength for oxygen can be calculated using 
the following empirical formulas (5): 

y = >o [Pd + 1.3e] °‘ 9 [GHz] (4-1.13) 

S = S 0 Pd (I-) 3 exp [d(l - ^) ] [KHz] (4-1.14) 

where Y 0 , S 0 , and d are constants for individual absorption lines (see 
Table 1). 
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For water vapor the line width and line strength have similar form- 
ulas: 

v- ro [Pd + 4.8ej °‘ 6 [GHz] (4-1.15) 

S = S 0 e exp [d (1 - ^)] [KHz] (4-1.16) 

B. Statistics 

To examine the variation of the atmospheric microwave refractivity 
at the earth's surface, the refractivity was calculated for each meteor- 
ological observation using the expressions developed in Section A and 
the probability density functions, diurnal variations, and seasonal vari 
ations were plotted. As was expected, the atmospheric refractivity var- 
ies diurnally and seasonally. 


1 . The Frequency Independent Refractivity 

An approximate equation which yield? values for N 0 within 0.02 per- 
cent of those obtained using Equation (4-1. 7 ) is (7): 


No 


77.6 

T 


( p + 


4810 


e) 


(4-2.1) 


where P is the total pressure, i.e. 

P = Pd + e [mbar] 

Since pressure does not vary significantly (Figures 8 and 9) and tempera 
ture varies within 253-313° [°k] (Figures 2 and 3), the frequency inde- 
pendent refractivity is highly dependent on water vapor pressure. Thus, 
the shape of the probability density function of N 0 is similar to that 
of water vapor pressure (Figures 28 and 29). 
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Even though the diurnal variation of water vapor pressure is very 
small (Figures 18 and 19), the normal temperature rise and pressure de- 
crease in the afternoon causes N 0 to decrease during that period (Fig- 
ures 30 and 31). 

In the summer, the average diurnal variation is about 13 
[ppm]. The average monthly variation is about 40 [ppm] during the 
winter and 80 [ppm] during the summer. The maximum and minimum of 
the frequency independent refractivity during 1949-1958 and 1967-1976, 
are 385 [ppm] and 275 [ppm], respectively (Figures 32 and 33). 

Figures 34 and 35 show the diurnal and seasonal variation patterns 
of temperature, pressure, water vapor content, and frequency independ- 
ent refractivity in a composite display for comparison purposes. 

2. Frequency Dependent Refractivity 

At low frequencies, the refractivity is simply the frequency inde- 
pendent value which is the summation of the low frequency tails of all 
of the absorption lines. As the operating frequency increases, the ef- 
fects of the individual absorption lines are observed; the refractivity 
goes through dispersion variations as the operating frequency passes 
through the resonance frequencies of the absorption lines. We plot the 
near, deviation of the refractivity from its freauency independ- 
ent value in Figure 36. The dispersion variations at 22, 60. 118, 183, 
325, and 380 GHz are shown. Of the water vapor absorption lines, the 
lines at 68 and 321 GHz are so weak that no significant dispersion is in- 
troduced near those frequencies. Note that the dispersion variations at 
the water vapor absorption lines are very strong during the humid sea- 
son. 
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Figure 30. Averaae diurnal variation of N 0 for the period 1949- 
1958/ 
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Fiqure 31. Average diurnal variation of N 0 for the period 1967- 
1976/ 
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c iqure 33. Daily maximum, average, r.iim'num, and extreme N 0 by 
months for tr„_ period 1967-1976. 
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Figure 34. Average diurnal variations of temperature, pressure, 
water vapor pressure, and N 0 for the period 1949-1958. 
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The Imaginary part of the refractivlty Is shown as a function of 
frequency In Figure 37. Since the oxygen molecule's contribution to the 
imaginary part of refractivlty does not depend on water vapor pressure 
(see Equation (4-1.14)), Its value In the vicinity of the oxygen absorp- 
tion lines (60 and 118.75 GHz) Is not affected much t; the change of the 
atmospheric state (see Figure 37). For clarity of the plot only January 
and July have been shown In these cases. 

Since communications engineers are more Interested In attenuation 
than in the refractivlty, the attenuation variation has also been plot- 
ted In Figure 38. This may be used to estimate the variation of attenu- 
ation for the clear sky. 

Because of current interest in propagation effects at twenty and 
thirty GHz, refractivlty statistics were derived for these frequencies 
and are presented in Figures 39 through 44. Although the most probable 
values of specific attenuation at these frequencies are 0.C5 and 0.03 
dB/km, respectively, we note that these values can exceed Q.2 and 0.1 
dB/km, respectively, during the humid sumner season. The excursion at 

20 GHz Is much larger than at 30 GHz due to the strona influence of the 
22 GHz water vapor absorption line. 
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Figure 37. Average imaginary part of refracti vity. 
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PROBABILITY DENSITY 



PERIOD s 1 /I /49-12/31/58 
FREQUENCY s 20.0 GHZ 
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Figure 39. Probability density function of N" at 20 GHz; £N"= 

1 .6 x 10-3 [ppn] . 
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Figure 40. Average diurnal variation of N' ; at 20 GHz. 
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Figure 43. Average diurnal variation of N" at 30 GHz. 
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CHAPTER V 


ATMOSPHERIC REFRACTION 

In a standard atmosphere, pressure and water vapor decrease rapidly 
with height, while temperature decreases slowly with height from the 
ground level to the troposphere. It is evident in such an atmosphere 
that the refractive index of the atmosphere decreases approximately in 
an exponential fashion with height; and, as a result, waves travelling 
in the atmosphere are bent downward under normal atmospheric conditions. 
The effects of this atmospheric refraction of radio waves are usually 
expressed in terms of the angle through which a radio ray turns as it 
passes through the atmosphere. To the radar engineer, '‘t is important 
to estimate the elevation angle error which is a measure of the differ- 
ence between the apparent elevation angle as indicated by radar and the 
true elevation angle. Usually the calculation of the effects of the at- 
mospheric refraction depends on the refractive index profile model used 
and the numerical integration techniques. 

In this chapter, we will develop analytic expressions for the bend- 
ing angle, elevation angle error, and range error in an exponentially 
tapered spherical atmosphere. 

A. Refractive Index Profile 

In atmospheric refraction studies, it is necessary to have a know- 
ledge of the behavior of the refractive index as a function of position. 
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It Is conrnon to assure that the refractive Index of the atmosphere Is 
spherically stratified with respect to the surface of the earth. Neg- 
lecting the effect of horizontal variations seems quite reasonable in 
the tropospheric region because of the relatively slow horizontal change 
of refractive index in contrast to the rapid decrease with height under 
normal conditions. 

An examination of long term climatic data of N $ and of AN, where 
aN is the decrease in refractivlty over the first kilometer in the at- 
mosphere and N $ is the refractlvity at the surface of the earth, reveals 
that a high correlation exists between these two parameters. The empir- 
ically derived relationship (7) is, 

AN = 7.32 exp [0.005577 N $ ] (5-1.1) 

This relationship between aN and N $ suggests that a model of the atmos- 
pheric refractivity structure in the first kilometer above the earth's 
surface may be obtained from N $ alone. A consideration of the standard 
atmosphere (Figure 1) and Equation (4-1.7) indicates that the refractiv- 
ity generally decreases exponentially with height. 

The exponential model specified by assuming a single exponential 
distribution of N is 

N = N $ exp {-c (r-r s )} (5-1.2) 

where N s = refractivity at radius r s 



Even though this exponential model does not givea good representation of 
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refractlvlty above 5km, It Is widely used In theoretical studies be- 
cause of the advantage of being an entire function. Furthermore, since 
most of the refraction effects occur In the lower portion of the tropo- 
sphere where the refractivlty Is largest, the accuracy of the model for 
larger heights is of less importance. 

If we assume « 1 (typically ^ = 0.13), then C = - In 

^1 - . Hence, our simplified exponential model will be 

described as: 

N = N s exp (5-1.3) 

This exponential model and Equatior. (5-1.1) will be used In the follow- 
ing section to develop the bending angle. 

B. Bending Angle 

If it is assumed that the refractive index is a function only of a 
height, then the path of a radio ray (Figure 45) will obey Snell’s law 
for polar coordinates: 


n^r^ cos G-j = nr cos e = n cos 


(5-2.1) 


The classical expression for the angular change, t, or the bending of a 
ray travelling frcm a point where the refractive index is n^ to a second 
point where the refractive index is n ? is given by (7) 



(5-2.2) 


where e is the local elevation angle. Tne differential equation for t 
will be: 
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Figure 45. Geometry of the refraction of radio waves. 
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d T = - cot G dn 
n 


(5-2.3) 


Since n ranges from about 1.0000 at very high altitude to about 1.0003 
at the surface of the earth. Equation (5-2.1) can be approximated as 
cos 0-j = r cos e 

cr 

r ! 

cos e = — cos 9, (5-2.4) 

r 1 

Now, by using trigonometry. 


cot e = 



(5-2.5) 


Substituting Equation (5-2.5) into Equation (5-2.3) yields, 
d T = - 


r^ cos O 


V 


2 2 2 
r - r^ cos &j 


1 . dn 

- -7— dr 

n dr 


(5-2.6) 


1 


and, since — ^ 1, 


dr x 


r^ cos 


V" 2 • r i cos! 


or 


dn 

dr 


dr 


T " 


r.| cos 0-j 


1 


V 


r^ cos ©^ 


dn 

dr 


dr 


(5-2.7) 


If the rate of change of the index of refraction with radius, 
is known, the bending angle, T , may be calculated. In Section A, we de- 
veloped an exponential model for refractive index structure: 
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(5-2.8) 


N - N $ exp [- (r-r s )] 

So, the derivative of refractive Index Is 
^ = - aH exp [- j|~ ( r - r s )l x 10 6 


where 

aN = 7.32 exp [0.005577 N s ] 

N s = refractivity at the surface of earth 
r s = radius of the earth 

Combining Equation (5-2.9) with Equation (5-2.8) gives 
t = r^ aN ' 10“ 6 cos e 1 exp [ r $ ] 


exp 


r - 


AN 



r ] 


dr 


-*.2 2 
r 2 - r i cos e-j 


Let us now change variables, h = r - r^, so that h represents t 
above the starting point, r^. And let H = r^ - r^ be the total 
in height when moving from point 1 to point 2. The Equation (5 
will be 


t = r^ aN • 10" 6 cos e 1 exp [- ~ (r^ - r^)] 



exp (- ^ h) 

-^/(h + r^ 2 - r 2 cos 2 0 1 


dh 


(5-2.9) 


(5-2.10) 

he height 
increase 
- 2 . 10 ) 


(5-2.11) 
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Now we note that, even for a large H, the Integrand is forced to zero by 
the exponential in the numerator. It is apparent that almost all of the 
refraction occurs within the lowest twenty km of the atmosphere. With 
h = 19km, the numerator, exp (~j^0» drops to 0.1. So the second order 

term in the denominator of the integrand may be ignored: 
t = AN * 10~ 6 cos 0-j exp [- jp- (r.j - r $ ) ] 


f H •» H7 h) 

J 0 sin 2 e, 


.. dh 


(5-2.12) 


h 2 2 

Introducing another change of variables, 2^— + sin c 0, = x , gives: 

r l 1 


i = r 1 an • 10' 6 cos 0 1 exp [- ^ (r ] - r $ )] exp [ j- sin 2 0-j] 


'2H . 2 

j=P- + sin e. 


exp [- 


AN r l 2- 

jqr x ■ 


(5-2.13) 


AN r l 2 _ J. 

» M X " * » 


And, with yet another change of variables, p- p- x 
(5-2.13) reduces to: 


Equation 


t * /’2rj AN ~ N^ 10" 6 cos 0j exp [- jp ( r p r s )l exp [jp j- sin 2 


IjT ^ H+ <T sip2 'M 

r sin_ 0 i~ 


exp (-t ) dt 


(5-2.14) 
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Since the error function Is defined as 
we have 


* <x) n JT 



exp (-x ) d x 


« 


= y I r ! AN • N s ' 10" 6 cos 0 j exp [- ^ (^1 -r s )] exp [|^ 


sin 2 ©i] 




(5-2.15) 


This Is the general expression for the bending angle in a single expo- 
nential model defined by Equation (5-1.3). 

We can calculate the bending angle for various situations by using 
the following rational approximation of the error function (8): 


• to’ 1 - [Ha,xV -!•- +C<X) <5 ' 2 ' ,6) 

where 

i e ( X ) ! 5 3 X icf 7 

a 1 = 0.0705230784 a 2 = 0.0422820123 

a 3 = 0.0092705272 a 4 = 0.0001520143 

a 5 = 0.0002765672 a fi = 0.0000430638 

But for a large initial elevation angle, ©i, the error functions con- 
verge to 1 and sin c ©i becomes a very large value. So in this 

case, we must use the asymptotic relation for the error function (8): 


$( x ) = 1 - 


•/n x 

The criteria for this will be: 


exp (-x ) 


4 


jf- 2^ S1n 01 >> 


(5-2.17) 
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This condition corresponds to 0 ! » 2.89° for 'leal values of aN, N $ 
and r,. Then, the bending angle for high inlt*-/ elevation angle Is glv 
en by: 


•t = N c • 10‘ 6 cos 0 j exp [- (rj - r )] 

s 



exp [- H] 


+ sin 2 ©! 


, 0i >> 3° 


(5-2.18) 


Figures 46 and 4 7 show the bending angle versus refract ivity for 70km 
a"d 35,785km heights. The nominal height of a geosynchronous satellite 
above the earth of radius r s = 6,378 km is 35,785 km (9). Since most 
of the refraction occurs in the region below twenty km, there is no no- 
ticeable difference between the?, two figures. 


C. Elevation Angle Error 

For earth- satell ite communication links, where narrow beamwidth an 
tennas are often used, the elevation angle error, e, cannot be ignored. 
To find this error, consider the geometry shown in Fiqure 48. Surrming 
all the angles about the quadrilateral defined by t\ j radii rj and r- 
and two local tangents at points (T) and (2)g1ves the relationship: 

$ = T + 02 “ ©1 (5-3.1 ) 

From the geometry shown in Figure 48, 


tan s = — ' 
©2 


cos e - r 2 
s Tn <t> 


(5-3.2) 


Since e = G] - e, substituting Equation (5-3.1) in Equation (5-3.2) 
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Figure 47. Bending angle versus ref racti vi +y for various Initial 
elevation angles; H »■ 35,785 km. 
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Figure 48. 



Geometry for elevation angle error. 
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yields: 

tan * = ~ tanB 

T + ^onejtane 

_ -r 2 cosxcos 0 2 + r 2 sintsine 2 ♦ r iCOSG 1 
r 2 sinxcose 2 + r 2 cosxs1ne 2 - r 1 sin8 1 

From Equation (5-2.1), 
r, 2 

‘'xCosgi * — r 2 cose 2 


(5-3.3) 


(5-3.4) 


Dividing denominator and nianerator of the right hand side of Equation 
(5-3.3) by rjcosei dives: 


tan 


e - 


cost + sinttane 2 + n 2 

n l 


.5* 


sint + cosTtaner 


•■anG! 


or 


2 

/-cos T + sinTtar.o? + — 


e = Arctan 


\^sin 


t + cosnane 2 - 


tan 


’1 


"l 


.1 

/ 

/ 


(5-3.5) 


Here, the bending angle, t, can be obtained from the expressions in Sec- 
tion A, and the angle between the local ray tangent and the local hori- 
zontal at point 0, £’ 2 » can be obtained using Equation (5-3.4) if we 
know r lf r 2 , and Oj . 

In Figures 45 and 50, w° plot the elevation angle error for 70 and 
35,785 [km], respectively. Comparison of these figures shows that the 
elevation angle error Increases with the heigh* of the end of the ray 
path. 
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Figure 50. Elevation angle error versus refractivity for various 
initial elevation angles: H = 35,785km. 
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D. Range Err o^ 

When a radio signal travels from point (]) to point (?) (see Figure 
48), the time delay is longer than would be expected if the signal trav- 
eled through a free space for tvo reasons. First, the path is longer 
due to the ray bending. And, second, the velocity of propagation is low- 
er in a refractive medium. The true slant distance, R true , ma l r be ob " 
tained from the geometry of Figure 48. 


! tru e * V r i 2 + r 2 2 ' 2r i r 2 COS$ 


(5-4.1) 


-1 

Since <p = x - e x + cos ( — cosQj), we can calculate Rtrue we know 

ri, r 2 » and e*. The bending angle, x, may be obtained from Equation 

(5-2.15) or (5-2.18). The apparent distance is given by: 

•r 


R a = cT = 


nds 


(5-4.2) 


w^ere T is the time delay for the signal travelling from point (j) to 


0 and c is the velocity of light in free space. Since ds = dr 

and n = 1 + N x 10”^, the above equation will be: 


-r 


nk dr * 10-6 


N 


sine 


dr 


(5-4.3) 


Thus, the total apparent range error, aR = R a - Rtrue’ 1S 9 iwen by: 

•r 

T Ir 

(5-*».4) 


aR 


^ 1 -6 

Jr x :W dr + 1 ° 


= / 


N 


J r l 


sine 


dr - R 


true 


The first term represents the length of the curved ra/ path. Since the 
difference in length between the straight path and the curved ray path 
does nnt represent a significant portion of the total range error except 
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at very small initial elevation angles, below about 3° (7), the range 
error for elevation angles greater than about 3° may be qiven by 


aR 2? 10“ 6 f 2 4-r 

J r i sin© 


dr 


(5-4.5) 


Substituting Equation (5-1.3) into Equation (5-4.5), and changing vari- 
ables, r - ri-»h and r : - r^+H, yields: 

aR = N. x lcf 6 exp [- (rj . r ,] 


(1+ £-) exp [- £ h] 
H r i N s 




dh 


(5-4.6) 


COS 2 ©! 


Note that the integrand is significant only when h is very small. Neg- 
lecting the higher order terms in , and from Equation (5-2.12), the 
range error can be related to the bending angle, x, by 
N. 


= - s - T 
A M 


(5-4.7) 


AN cose; 

The range errors versus the surface refractivity are plotted in Figures 
51 and 52 for H = 70 km, H = 35,785 km, respectively. Tnese may be used 
for the first c-der approximation of the range error. It should be em- 
phasized that the difference between the straiqht and the curved path 
contributes nearly thirty percent of the total range error at zero ele- 
vation angle (7). 
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Figure H. Range error versus reflectivity for various initial ele- 
vatio.i angles; H= 70km. 
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CHAP't eft VI 
SUMMARY 

The atmosphere is a most changeable propagation media; and its 
electrical characteristic, microwave refractivity, is i>oendent on temp- 
erature, pressure, and water vapor content. The statistics of these 
meteoroloaicai quantities at Columbus, Ohio, were generated using two 
ten-yea»' period weather data bases from the National Weather Service. 

The microwave refractivity was calculated using these meteorologi- 
cal data bases and Its variations were plotted and analyzed. In the 
calculation of the microwave refractivity for the EHF range, the G'-oss 
line shape was used. The frequency independent refractivity varies from 
about 280 to 380 [ppm] and the average seasonal variation is about 35 
[ppm]. The average seasonal variation of attenuation due to the atmos- 
pheric gases is about 0.13 dB/km at 20 GHz, and 0.0b dB/km at 30 GHz. 

The maximum variations cf these are 0.2 and 0.1 db/km, respecti vely. 

3ased on a single exponential model for the atmospheric refractiv- 
ity structure, bending angle, elevation angle error, and range error 
were developed and plotted. Most of refraction occurs in the lower re- 
gion of troposphere because the refractivity decays exponentially with 
height. For small elevation angles, below about 5°, the refraction ef- 
fects are very sensitive to the refractivity at the earth's surface. 
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APPENDIX A 


CODING AND FORMAT OF METEOROLOGICAL DATA 


Character Set 


HEXADECIMAL OCTAL 


1 

FI 

01 

2 

F2 

02 

3 

F3 

03 

4 

F4 

04 

5 

F5 

05 

6 

F6 

06 

7 

F7 

07 

8 

F8 

10 

9 

F9 

11 

0 

FO 

12 

a (blank) 

40 

20 

- 

60 

40 

* 

5C 

54 

& 

50 

60 

( A 

Cl 

61 

B 

C2 

62 

C 

C3 

63 

D 

C4 

64 

< E 

C5 

65 

F 

C6 

66 

G 

C7 

67 

H 

C8 

70 

. i 

C9 

71 

fj 

D1 

41 

K 

D2 

42 

L 

D3 

43 

1 M 

D4 

44 

S N 

D5 

45 

0 

D6 

46 

P 

D7 

47 

Q 

08 

50 

ft 

D9 

51 

* 

CO 

72 


DO 

52 


EO 

32 
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* CO through C9 in hexadecimal represent positive decimal numbers. 

** DO through D9 in hexadecimal represent negative decimal numbers. 


Record Format 





TAPE 

TAPE 


FIELD NUMBER 

POSITIONS 

ELEMENT 

001 

001 


004 

Tape deck number 

002 

005 

- 

009 

Station number 

003 

010 

- 

on 

Year 

004 

012 

- 

013 

Month 

005 

014 

- 

015 

Day 

101 

016 

_ 

017 

Hour 

102 

018 

- 

021 

Ceiling height and indicator 

103 

022 

- 

025 

Horizontal visibility and indicator 

104 

026 

- 

027 

Wind direction - 16 points 

105 

028 

- 

030 

Wind speed 

106 

031 

- 

033 

Dry bulb (air) temperature 

107 

034 

- 

036 

Wet bulb temperature 

108 

037 

- 

039 

Dew point temperature 

109 

040 

- 

043 

Relative humidity and indicator 

no 

044 

- 

048 

Sea level pressure 

111 

049 

- 

052 

Station pressure 

112 

053 

- 

057 

Sky condition and indicator 

113 

058 



To^al sky cover 

114 

059 



Total opaque sky cover 

115 

060 



Amount of lowest cloud layer 

116 

061 



Type of lowest cloud or obscuring pheno- 





mena 

117 

062 

- 

064 

Height of base of lowest cloud layer or 





obscuring phenomena 

118 

065 



Amount of second cloud layer 

119 

066 



Type of cloud - second layer 

120 

067 

- 

069 

Height of base of second cloud layer 

121 

070 



Summation amount of first two cloud lay- 





ers 

122 

071 



Amount of thi . d cloud layer 

123 

072 



Type of cloud - third layer 

124 

073 

- 

075 

Height of base of third cloud layer 

125 

076 



Summation amount of first three cloud 





layers 

126 

077 



Amount of fourth cloud layer 

127 

078 



Type of cloud - fourth layer 

128 

079 

- 

031 

Height or' base of fourth cloud layer 

129 

032 



Occurrence of thunderi orm, tornado, or 





squal 1 

1 30 

C83 



Occurrence of rain, rain showers, or 


freezing rain 


89 



TAPE 

TAPE 


FiaO NUWER 

POSITIONS 

tLEPCNT 

131 

084 

Occurrence of rain squalls, drizzle, or 
freezing drizzle 

132 

085 

Occurrence of snow, snow pellets, or ice 
crystals 

133 

086 

Occurrence of snow showers, snow squalls 
or snow grains 

134 

087 

Occurrence of sleet, sleet showers or 
hail 

135 

088 

Occurrence of fog, blowing dust, or 
blowing sand 

136 

089 

Occurrence of smoke, haze, smoke and 
haze, dust, blowing snow, blowing 
spray 

137 

090 - 091 

Wind direction - 36 points 

138 

092 - 094 

Blank 

139 

095 

Record mark 


201 - 239 096 - 175 Second observation \ These observations 

30! - 339 176 - 255 Third observation follow the same 

401 - 439 256 - 335 Fourth observation » format as fields 

501 - 539 336-415 Fifth observation 101-139 (Tape Posi- 

601 - 639 416 - 495 Sixth observation tions 016-095). 
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TAPE 

TAPE 


TAPE 

CODE DEFINITIONS 

FIELD NUMBER 

POSITIONS 

ELEMENT 

CONFIGURATION 

AND REMARKS 

001 

001 - 004 

Tape deck number 

1400 - 1499 

Used to distinguish different data 
sources. See current list at beginning 
of manual . 

002 

005 - 009 

Station number 

01001 - 98999 

Unique number used to Identify each 
station. Usually a WBAN number but oc- 
casionally a WHO number. 

003 

oio - on 

Year 

00 - 99 

Year of observation. 00-99 * 1900-1999 

004 

on - 013 

Month 

01 - 12 

Month of observation. 01-12 * Jan. -Dec 

005 

014 - 015 

Day 

01 - 21 

Day of month 

101 

016 - 017 

Hour 

00 - 23 

Hour of observation In local standard 


time. 00-23 « 0000-2300 LET 


106 


109 


111 


031 - 033 Dry bulb (AIP) tem- 
perature 


041 - 043 Relative humidity 
« 

049 - 052 Station pressure 


OOl - 130 Specified temperature In whole degrees 
fahrenhelt. 

001-130 - -1° - -1 30°F 

00$-14$ * 0° - +140°F 

a aa * Unknown 

aa* * Original value Invalid 

001 - 100 Relative humidity In whole percentage 
aaa, A a* aaa * Unknown 

aa* * Original value invalid 


1900 - 3999 Pressure at station level In Inches 
a aaa, aaa* and hundredths of hg. 

1900-3999 = 19.00 - 39.99 In hg. 
aaaa * Uknown 
aaa* - Invalid 



